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EXECUTIVE SUMMARY 


INTRODUCTION 

Earth and earth-space communications links are affected by 
precipitation along the transmission path. Terminals set up to take 
advantage of space diversity must be spaced in a manner dependent on 
rain intensities, rain cell sizes, and spaclngs. Unwanted Interference 
also occurs because of scattering from the precipitation particles. 
Communications engineers urgently need statistical information on rain 
shower intensities and their geometrical properties. 

A NASA sponsored experiment was carried out at Wallops Island, 
Virginia in 1973 to provide such data. One of the S-band radars, the 
Space Range Radar (SPANDAR) , was modified to meet the needs of a radar 
specifically aimed at studying precipitation processes. 

This report. Volume I, describes the radar changes, the rain 
cell statistics program, and the data collected. In Volume II are found 
details of the data reduction and analysis program. Volume III contains 
some of the more Important memoranda Issued during the course of the 
experiment . 

SYSTEM DESCRIPTION . 

Radar 


SPANDAR is a high-power, narrow-beam radar originally designed 
for long-range tracking.. Fine resolution is achieved by the 60 ft an- 
tenna (0.4 deg) and the 1 ysec pulse length. Great care was taken with 
pre-experiment radar measurements and calibrations during the experiment 
to reduce the combined radar and processing errors to less than 2 dB. 

There were changes in the original radar configuration occa- 
sioned by the specific needs of the rain cell experiment. These modifi- 
cations were: frequency diversity capability, digital integration, 

computer control of the antenna mount, extended linear range of the 
receiver, and digital and analog video recording providing on-line as 
well as off-line data processing capability. As a consequence of these 
modifications the radar can be considered to approach the ideal radar 
for the radar meteorologist. 

Frequency Diversity 

It is well known that since meteorological targets fluctuate 
from pulse to pulse it is necessary to average a number of statistically 
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Independent pulses to achieve a reliable measure of reflectivity. The 
largest number of pulses averaged is generally limited by scanning con- 
siderations while the shortest meaningful time between pulses to be 
averaged depends on their tlme-to-independence . For typical weather 
targets at S-band the tlme-to-independence is about 20 msec. 

Shifting the r-f frequency from pulse to pulse is an excellent 
way to obtain independence in a short time. With the 1 psec pulse it 
was necessary to jump in frequency at least 1 MHz for say 20 successive 
pulses before returning to the original frequency. The actual frequency 
shift was in an interlace mode where the frequency interval between 
pulses was 12 MHz. The interlace pattern was selected to minimize 
second- time-around weather targets; such interference can be a problem 
when the rain is heavy. For a 12 MHz frequency jump the second-tlme- 
around targets fall well outside the receiver bandwidth. 

To accomplish the goal of providing statistical information on 
rain cells one needs a three-dimensional snapshot of the cell structure. 
This can be achieved by making a raster scan through the cell before the 
cell changes. A scan rate of 3 deg per sec was considered sufficiently 
fast to do this. This rate is adequate if one requires one sample 
reflectivity per beamwldth and every successive pulse is Independent of 
its predecessor. It can be seen that the appropriate scan rate is a 
compromise between scanning fast enough to see 'the entire cell structure 
before it changes and slow enough to ensure accuracy of measurement; 
frequency diversity allowed sampling and averaging 128 pulses for 1 
sample in each beamwldth. 

Digital Integrator 

An 871 bln digital Integrator capable of integrating 128 to 
1024 video pulses was added to the system. It is true that for a sta- 
tionary process one obtains a more accurate measure of reflectivity the 
larger the number of pulses in the average. However, adding more sam- 
ples after about 50 yields an increase in accuracy only very slowly. 

The specific number of 128 pulses in the average was arrived at for prac- 
tical reasons . 

Computer Control of Antenna 

In contrast with most of the previous research performed with 
the Wallops radars, where manual operation was adequate, this experiment 
required that antenna scan sequences be executed with little loss of 
time between sweeps. A PDP-8 computer was Interfaced with the antenna 
mount to perform sequences rapidly and automatically. Stored programs 
for various scan sequences can be called up quickly via teletype control 
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resulting in more-or-less stationary rain cell snapshots. Complete 
three-dimensional snapshots of cell structures were scanned within a 
four-minute interval; this included azimuthal sectors up to 60 deg in 
width. 

Extended Dynamic Range Log Amplifiers 

For this experiment it was deemed important to measure weak 
rain areas as well as the heavy showers. Calculations for SPANDAR 
showed that in order to measure a 1 mm per hr rain at the longest ranges, 
140 km, as well as heavy thunderstorms at close ranges, 10 km, it would 
require a receiver linear over a 90 dB dynamic range. This assumed 10 
dB pulse-to-pulse fluctuations of the received weather signal. A new 
solid state log IF amplifier was designed and added to the system to 
extend its linearity over this range. 

Digital and Analog Video Recording 

An on-line recording system was made part of SPANDAR for the 
rain cell experiment. The output of the video Integrator was digitized 
and recorded on a 9-channel, 1600 bits per in. synchronous tape recorder 
operating at 45 in. per sec. In addition to the Integrated output from 
each resolution cell, the digital tape stored time, antenna mount po- 
sition, processing parameters as well as all pertinent radar parameters 
including pulse-to-pulse calibration signals. 

An off-line system also was Incorporated into SPANDAR as a 
precaution against loss of data. The basis of this system was an Ampex 
FR-910 video tape recorder. It is a four-head, spiral-scan recorder to 
record radar video on a frequency-modulated subcarrier with an input 
signal bandwidth of 6.5 MHz. Two audio channels are available for re- 
cording voice annotation and/or multiplexed data ancillary to the video 
signals . 


Radar video was recorded on the FR-910 and Integrated video on 
the digital tapes. If any malfunction occurred in the processor the 
original radar video was replayed later through the processor to simu- 
late real-time recording. In fact, the original video tapes still exist 
and can be replayed at any time for other purposes. The video tapes can 
also be played back through RHI and FPI scopes using the ancillary data 
by using the special interface equipment designed for this purpose. 

EXPERIMENT DESCRIPTION 

The experiment was carried out in two parts. The first part, 
during the fall of 1972, was essentially a check-out period for the 
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equipment and procedures, but some useful data were obtained. During 
the summer of 1973 the full-scale experiment took place using all the 
equipment as described here. 

The general purpose of the experiment was to determine the 
statistics of intense rain cells using radar. In particular, the ob- 
jectives were to (1) determine the statistics of the field of rain cells 
in terms of such descriptors as cell size, intensity, spacing, number 
density, and altitude dependence; and (2) find the detailed structure of 
individual cells leading to a model of an individual cell for use by 
communications scientists and meteorologists. 

For the first objective the operational mode was to take a set 
of PPI scans, usually over 360 deg, at a 2 deg elevation angle at one- . 
half hour intervals. For the cell model objective a mode of operation 
was used involving rapidly-taken PPI raster scans. In this mode a 60 
deg sector was selected and one, or a series of, raster scans were made 
to examine the three-dimensional structure and its changes with time. 

The preliminary experiment ran from 15 September to 6 October 

1972. The full-scale experiment covered the period 9 May to 31 August 

1973. 

METEOROLOGICAL CONSIDERATIONS 

As is apparent in any radar-weather experiment certain synop- 
tic mesoscale and perhaps microscale weather data are needed. One must 
forecast precipitation in order to carry out the experiment economically. 
One needs point rainfall measurements and other meteorological observa- 
tions to help relate the radar measurements to the broader meteorological 
situation and also to help in extrapolating the rain cell statistics to 
other climatic regions. 

The Weather Service at Wallops Island was used as the primary 
source of weather information. Surface observations and radiosonde data 
were made available along with all the teletype and facsimile charts 
from the networks. In addition, weather briefings were held each morning 
which Included local forecasts. A tipping bucket rain gage was installed 
about 12 mi north of SPANDAR especially for this experiment; this was 
augmented by rainfall data at Wallops, Patuxent, Dulles, Atlantic City, 
and Richmond. 

ACCURACY AND CALIBRATION 

Considerable effort was placed on measuring the radar system 
carefully. The antenna gain is now known precisely, system stability 
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was checked by frequent calibrations. The overall SPANDAR accuracy may 
be stated to be within 1.7 dB standard deviation. 

DATA SUMMARY 

In the 1972 experiment there were 10 rain days and 123, 360 
deg PPI scans at one-half hour intervals and 35 PPI raster scans were 
obtained to study individual cells in three dimensions. In the 1973 
experiment there were 33 rain days with 96 one-half hourly PPI scans and 
403 PPI raster scans. A total of 1141 discrete rain cells were processed 
from these data. 

Fortunately for this experiment, the 1973 summer period had 
more rain and more days of rain than the preceding 10-yr average pri- 
marily because May and August exceeded the averages by a substantial 
amount . 
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INTRODUCTION 


Ground stations transmitting to communications satellites will 
under certain propagation conditions suffer atmospheric attenuation and 
may also Interfere with beyond-the-horlzon ground communications links 
operating on the same frequencies. Because of severe spectrum limita- 
tions and a consequent Intense Interest In frequencies above 10 GHz, It 
has become necessary to obtain propagation and Interference Information 
on which one can base decisions In communications systems design. One 
must either limit transmitted power or assure adequate separation of the 
stations. In either case, decisions must be based on quantitative 
Interference measurements or accurate predictions based on such measure- 
ments . 


There are various mechanisms for propagating short electro- 
magnetic waves beyond the horizon. For the radio frequencies being dis- 
cussed here, the Important causes of scatter are ducts, layers charac- 
terized by high variance In refractive Index, and discrete scatters such 
as aircraft, birds, and precipitation. Crane (Ref. 1) has concluded 
that precipitation scatter Is the most Important of these mechanisms for 
the communication Interference problem. 

In that same paper he also concluded that the statistics of 
the small, strongly reflecting rain cells were unavailable. Such cells 
are difficult to describe statistically because of their Infrequent oc- 
currence. Ground-based measurements obtained with rain gages are In- 
adequate for this purpose because of the small likelihood of a cell 
passing over a gage. Only a large network of closely spaced gages can 
be used to study such cells, for example Semplak and Keller (Ref. 2). 

The most appropriate tool for observing such phenomena Is high 
resolution radar. Large areas (~4 x 10** km^) can be observed frequently 
(within minutes) and with high resolution (~.4 km). Statistics such as 
the distribution of cell sizes and Intensities, typical cell separation, 
frequency of occurrence and cell lifetimes are but a few of the Important 
characteristics measurable with radar. A reasonable approach Is to 
determine the above statistics using radar and relate them to local rain 
gage measurements and a characterization of the rain type. It Is likely 
that these relationships can ultimately be Inverted to predict the 
expected statistics of rain In other areas using only local surface 
rainfall data. 

The objective of this experiment was to collect radar data 
which could yield such statistics on the spatial distribution of summer 
rain cells. Ideally one would like simultaneous, fine-scale measure- 
ments In both time and space over a large volume, l.e., over a 200 km 
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radius and all altitudes to 15 km. In practice, however, this Ideal set 
of data cannot be realized using only a single radar for a variety of 
practical reasons. For example, because of the long tlme-to-lndependence 
of meteorological scatterers and the number of Independent radar samples 
on target needed to provide a good statistical estimate of reflectivity, . 
the time to probe such a large volume with a 0.5° beam Is much longer s 
than the fine time scale needed to establish cell lifetimes. This r 

incompatibility leads to the definition of two well defined objectives, 
each of‘ which Is capable of being handled by a radar. These objectives, - 
then, defined particular modes of operation for the radar system in the 
presently described experiment. 


■ The first objective concerned the determination of the sta- 
tistics of the field of rain cells, i.e., sizes, and spacings, the so-, 
called 'cell field mode. The second objective was to determine the 
detailed structure of individual rain cells leading to the description , 
of a cell model. In this case, relationships such as reflectivity as a 
function of altitude, size (both vertical and horizontal) and time 
development were needed. This mode of operation involved probing se- 
lected rain cells at a fine time scale and with fine spatial resolution. 


The experiment was carried out in two phases. The first 
occurred during September and October 1972 and emphasized testing and 
calibrating the equipment. The second phase was the full-scale data 
collection period providing the statistical base of one summer season, 

May through August 1973. 

Background for this experiment is contained in Refs. 3 through 
8. These technical notes contain the details of the equipment, the ex- 
perimental procedures and the original data analysis plan for obtaining 
the statistics of the time-and- space distribution of rain. 

This final report is divided into three parts. Volume I con- 
tains a description of the radar and associated equipment and the modi- 
fications made to carry out the experiment, certain meteorological con-: . 
sideratlons and a summary of the data collected during the experiment, 
The^data reduction, analysis techniques and results are presented in 
Volume'll. Some of the more important Technical Notes are reproduced - 
in Volume III. Also, since this experiment was so multi-faceted, it v 
was deemed advisable to prepare a cross-index to aid in finding source 
material in the Technical Notes; this cross-index wlll.be found in 
Section 8 of this volume. , ,,, 
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SYSTEM DESCRIPTION 


2.1 SPANDAR 

The Applied Physics Laboratory of The Johns Hopkins University 
has been operating a Radar Atmospheric Research Facility (RARF) at 
Wallops Island continuously since 1965. Support for the program came 
jointly from the Air Force and NASA. The radar facility is owned by 
NASA, Wallops Island and is part of their comprehensive radar, rocket 
and weather facilities. To date, RARF has been used in a variety of 
research problems associated with clear-air atmospheric processes as 
well as with precipitation. 

The Spa^e Range Radar (SPANDAR) is the S-band radar system 
used for the rain statistics problem being discussed here. It is a pre- 
cision, long-range, conical scan, tracking radar used primarily for the 
collection of deep space trajectory data and secondarily for range 
safety, orbit determination, and vehicle performance analysis. Its 
special features include a pulse doppler tracking subsystem and a digi- 
tal range subsystem that yields unambiguous range data to 250,000 nmi. 

The SPANDAR facility is located on the mainland, about 2 mi 
west of Wallops Island. The facility building houses transmitters, 
receiver, ranging equipment, data handling equipment, recording equip- 
ment, and a complete test and repair shop. Adjacent to the building is 
the antenna tower housing the receiver, rf-head equipment, parametric 
amplifiers, encoders, antenna drive components, and antenna pedestal. 

The technical characteristics of SPANDAR are shown in Table 2.1. 
Note that two transmitter tubes are available. 

One transmitter tube was designed by Raytheon and is unique. 
Replacement parts must be built as needed so that cost of repair and 
replacement is fairly expensive. The other is from an FPS-18 radar. 

FPS-18 radars are used extensively so that replacement parts are avail- 
able in quantity at low cost. In addition, maintenance procedures are 
well established. Peak power output of the FPS-18 transmitter is 1.2 MW 
as opposed to the 5 MW from the Raytheon tube. Since high power is not ■ 
necessary for the purpose of studying rain, the FPS-18 was considered 
adequate and was used throughout the experiment described here. 

The receiver system includes a parametric amplifier with ap- 
proximately a 20 dB gain. Saturation with the paramp operating is 
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Table 2.1 


SPANDAR Configuration Characteristics 


Transmitter Subsystem 


Frequency range 

2,700 to 2,900 MHz 

Peak power output 

Raytheon Transmitter 5 MW 
FPS-18 Transmitter 1.2 MW 

Pulse width 

1, 2 or 5 ysec 

Pulse repetition 

- - 

rate (PRF) 

160, 320 or 640 pps 

Receiver Subsystem 


Frequency range 

2,700 to 2,900 MHz 

Type 

Superheterodyne 

Dynamic range 

70 dB (modified to 90 dB) 

Noise figure 

2.5 dB 

Sensitivity 

-117 dBm for 240 KHz bandwidth 

I.F. frequency 

30 MHz 

Low noise device 

Parametric amplifier 

Bandwidth 

Matched to pulse width 

Antenna Subsystem 


Reflector 

60 ft parabola 

Feed 

Nutating 

Polarization modes 

Circular, vertical or horlzpntal 

Scan 

Conical 

Gain 

52.8 dB 

Beamwldth 

0.39° (half power points) 

Beam cross over 

1.5 dB 

Azimuth coverage 

360® continuous 

Elevation coverage 

0 to 90® 

Tracking rates 

5® sec“^ 

Angle slewing rates 

12° sec”^ azimuth; 17° sec“^ 
elevation 




around -50 dBm so that when the paramp is removed from the line, a 
saturation level of -30 dBm is reached. 

The mixer-preamplifier is located in the building at SPANDAR 
and the output of the paramp, which is located on the antenna feed, is 
brought down by cable. Higher level signals may be observed when the 
paramp is out by placing attenuators in the line prior to the mixer. 

The dynamic range of the log IF system was 70 dB, modified to 90 dB for 
this experiment, with a minimum detectable signal of -117 dBm with the 
paramp in the system. The SPANDAR antenna is a 60 ft diameter parab- 
oloid having a pencil beam whose width is 0.4 deg; its gain is 52.8 dB. 

A photograph of the antenna reflector is shown in Figure 2.1. It is cap- 
able of continuous azimuthal rotation at angular speeds up to 12 deg 
sec“^ in azimuth and 17 deg sec~^ in elevation. In elevation the antenna 
provides coverage up to 90 deg. 

There were changes in the original SPANDAR configuration occa- 
sioned by the specific needs of the spatial statistics experiment. 

These modifications were: 

a. Pulse-to-pulse frequency diversity into the FPS-18 
transmitter (Ref. 10) 

b. Computer control of the 60 ft antenna pedestal for 
various scan schemes 

c. Log amplifiers added to extend the dynamic range of 
the receiver 

« 

d. Integration of 128 pulses with digitization into 
871 range bins 

e. Computer control of the data acquisition and for- 
matting for processing, detection and analysis on 
an IBM 360/95 computer 

f. Digital and analog data recording on magnetic tape 

g. Pulse-to-pulse calibration of the radar video 
output 

h. Analog magnetic tape playback through the radar 
receiver and digital data acquisition system. 
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Figure 2.1. SPANDAR Dish and Feed 


Details of these modifications may be found in Refs. 9 , 10, 11, 12, 13 
and 14. 


2.2 FREQUENCY DIVERSITY 

Radar returns from meteorological targets fluctuate from pulse 
to pulse so that it is necessary to average many pulses to achieve a 
reliable measure of reflectivity. The pulses to be averaged must be 
statistically independent. A more detailed discussion of this point is 
contained in Section 3.3 on measurement accuracy. The greatest number 
of pulses averaged is generally limited by scanning considerations 
whereas the shortest meaningful time between pulses to be averaged de- 
pends on the time-to-independence, i.e., the decorrelation time. The 
time to independence, Xg, for typical meteorological scatters observed 
at S-band is about 10 to 30 ms (Atlas, Ref. 15 and Marshall and Hitsch- 
feld. Ref. 16). Samples taken at intervals less than Xg do not con- 
tribute significantly to the accuracy of the result although some bene- 
fit is gained by averaging partially correlated samples. The turbulence 
associated with thunderstorms would result in Xg being near the lower 
end of the above range. 

One practical method for achieving independent samples more 
quickly than x^ is by the use of frequency diversity (FD) . Wallace 
(Ref. 17) shows that independence can be achieved by pulse-to-pulse 
frequency shifting. The correlation coefficient p between two echoes Av 
apart in frequency is given by: 

2 

p _ r sin-n-xAv l 
L ttxAv j 

where x is the transmitted pulse width. Thus, complete decorrelation or 
independence from one pulse to the next is achieved if each pulse is 
shifted in frequency by 1/x. 

If no decorrelation technique is used such as frequency, spa- 
tial, or aspect diversity, an Inordinately long time is required to 
accumulate enough Independent pulses to calculate an accurate mean radar 
reflectivity. Frequency diversity (FD) capability was added to the 
SPANDAR radar for this purpose. A 1 ys pulse length was used so that 
the frequency had to be shifted by at least 1 MHz for 20 successive 
pulses or more before returning to the same frequency. 

The FPS-18 transmitter, used in this program, has a usable 
bandwidth of 30 MHz. There are two frequency stepping modes. In the 
first mode the frequency is stepped in 1 MHz Increments as described 
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above for 24 steps. The second is an interlace mode where the frequency 
interval between pulses is 12-13 MHz. This second logic, Figure 2.2, was 
included to minimize the second-time-around problem, i.e., the detection 
of targets wlt|iin the second (and even third) ambiguous range interval. 
The latter logic was used extensively through the 1973 experiment. 

The FD equipment consists of a solid state (tunnel diode) 
microwave oscillator unit controlled by digital frequency stepping 
logic, see Appendix A, Volume III or Ref. 14. The output of this oscil- 
lator is up-converted to the desired transmit and local oscillator 
frequencies by modulating the basic oscillator frequency (fo) with lower 
frequencies (240 MHz and 270 MHz) which are separated by the radar IF. 
Additional circuits included in the FD chassis to provide a calibrated 
reference signal are a time-delayed and pulsed sample of the transmit 
signal and a pulsed noise source. These two signals are combined in a 
microwave network and made available to the radar receiver in a common 
line. A block diagram of the FD unit is shown in Figure 2.3. 

The pulse delay and frequency-stepping logic generator unit is 
comprised of delay multivibrators, presettable counters, an adder and 
gates to provide the delayed pulses to the pulse modulator and noise 
diode source and frequency step circuit. The frequency step circuit 
accepts the delayed input trigger and steps the output frequency 8, 16, 

24 or 32 times before resetting to the starting frequency. A mode 
control switch is provided to select fixed frequency, linear steps 
(ramp) or an Interlace mode where an adder is enabled and disabled on 
alternate pulses to separate adjacent transmitted frequencies by the 12- 
13 MHz discussed previously. The digltal-to-analog converter and the 
shaping circuit convert the output of the digital frequency selector 
unit to an analog voltage, adjusts the gain for 1, 2, or 4 MHz per step 
and shapes the output voltage so that each digital step results in the 
same oscillator frequency shift. 

2.3 RADAR VIDEO PROCESSOR AND CONTROLLER 

Because of the requirement to obtain three-dimensional "snap- 
shots" of storm cells it was necessary to perform many of the data gath- 
ering and processing functions automatically. A PDP-8I computer was 
Installed to direct the antenna for the rapid scan sequences. Stored 
programs allowed the operator to select particular scan sequences for 
rapid execution. 

A digital integrator, also under the direction of the PDP-8I, 
integrated 128 video samples in each of the 871 range gates and recorded 
these data on IBM-compatible tapes. A photograph of the processing and 
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Figure 2.3. Block Diagram of Frequency Diversity Unit 


10 




















controlling equipment, the tape recorders and the teletype for input 
information is shown in Figure 2.4. 

A description of the control and processing equipment in terms 
of the method of processing on-line and off-line follows. 

2.3.1 On-Line Processing 

The on-line system performs four distinct functions; it digi- 
tizes and integrates the radar signal in each of the spatial resolution 
cells, records the integrated video on digital magnetic tape and provides 
antenna scan commands in real-time. The elements of this system shown 
in Figure 2.5 include the radar video processor (RVP) , digital tape 
recorders, the PDP-8I computer, the radar control and the RVP-Control 
Interface units. Refs. 18, 19, 20, 21. 

The function of the RVP is to sample the incoming radar signal 
in each spatial resolution cell, digitize these data and integrate for 
the required nvimber of pulses. To provide a reasonably accurate esti- 
mate of the reflectivity, the signal in each resolution cell must be 
averaged. A minimum of 128 pulses in the average was chosen which gives 
a statistical error of less than 0.4 dB, 80 percent of the time. The 
number of pulses Integrated was chosen based on considerations of prf, 
antenna scan rate and accuracy of the reflectivity estimate. A prf of 
960 sec”^ and 128 pulses in the average was used and the antenna scan 
rate was 3 deg sec“^. 

The spatial resolution in the radial (range) direction is one 
pulse length (1 ys which corresponds to 150 m) . The RVP accommodated 
871 range cells to cover the desired range increment from 10 to 140 km. 
The output for each cell is the mean log power received (log A^) in the 
form of an 8 bit word which is recorded on the digital tape recorder in 
IBM compatible format. Major features of the RVP are shown in Table 2.2. 

The digital tapes are 9 channel, 1600 bpi, synchronous machines 
which operate at 45 ips. Two tape transports are used with a controller 
for continuous operation. No time or data are lost due to changing tape 
reels. In addition to the integrated output from each resolution cell, 
the digital tape also stores time, the mount positions, processing 
parameters as well as all pertinent radar parameters including cal- 
ibration signals. Timing signals are provided to the RVP by the RVP- 
Control Interface. This interface also conditions the data from the 
radar to the RVP. The control of the RVP and digital tapes is provided 
by the PDP-8I computer. 
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Figure 2.4. Radar Video Processing Equipment Showing Digital Tape 
Recorders and Teletype Input 




Figure 2.5. Block Diagram of Radar Control Processing and Recording 
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Table 2.2 

Major Features of the Radar Video Processor 


RADAR VIDEO PROCESSOR 

871 Range Cells - 10 to 140 km 
Digitize Video, Integrate 128 Pulses 
Antenna Scan Rate 3 deg sec“^ 

Sample Video at 4 MHz Rate 
Output: 8-bit Word on Magnetic Tape 

MAGNETIC-DIGITAL TAPE 

9-Channel, 1600 BPI, 45 in. per sec 
Synchronous, IBM-Compatible 
Two Tape Transports for Continuous 
Operation 

Stores: Video, Time, Azimuth, and Elevation 

Processing Parameters 
Radar Parameters 
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The PDP-8I computer also controls the antenna mount through 
the Radar Control and the RVP-Control Interface. Programs are stored in 
the computer which command the mount to execute prescribed scanning 
sequences such as 360 deg scans. Range Height Indicators (RHI) or Plan 
Position Indicators (PPI) , raster scans, etc. Selection of scan param- 
eters such as azimuth start, stop, range segment to be recorded, execute 
commands, etc., are operator input to the computer via teletype. The 
computer also contains internal self checks which ensure proper operation. 

2.3.2 Off-Line Processing 

Whenever any element in the on-line data acquisition system 
was inoperable (the digital tapes and transports and that portion of 
the RVP-Control Interface which controls the RVP) , the off-line system 
was used. This system consists of an Ampex FR-910 video tape recorder 
and two interface units. Refs. 22 and 23. The tape recorder is a studio- 
type video recorder modified for instrumentation use. It is a four- 
head, spiral-scan recorder which records the video radar signals on a 
frequency-modulated subcarrier with an input signal bandwidth of 6.5 
MHz. Special amplitude and phase-balancing controls are provided to 
allow precise alignment of the signals to the four heads and eliminate 
switching transients. Two audio channels are also available for record- 
ing voice annotation and/or multiplexed data ancillary to the video 
signals. 


The record interface provided buffering of the incoming video 
signals to match the recorder's + 1.0 v input voltage and to mix the prf 
sync trigger on this video channel. Coded time (NASA format), bearing 
and elevation shaft encoder signals and other pertinent header (digi- 
switch entry) information was recorded as a multiplexed pulse train on 
one audio channel. The 960 prf is counted down to provide the shaft 
encoder strobe signal and the frame sync gate for each frame of the 
multiplexed input data. 

The playback interface had two functions. It allowed data to 
be played into an RHI or a PPI scope for inspection and it provided 
proper inputs for the RVP. Thus the interface unit had the capability 
of separating the prf sync pulse and the video signals and to restore 
the video to the original input levels. The transfer gain through the 
recorder was kept at unity. Major features of these two interface units 
are shown in Table 2.3. 

The digital data on the audio track is decoded for playback 
and analysis with selectable front panel light display, sin/cos outputs 
operate scope displays and a computer function organizes all the data 
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Table 2.3 


INTERFACE 


INTERFACE 


Major Features of Interface Units 

- RECORD 

Buffer to Match FR 910+ 1.0 v Input 

Mix PRF Sync Trigger on Video Channel 

Coded Time (NASA Format) ^ On Audio 

Azimuth & Elevation Shaft / Channel 
Encoder Signals j 

- PLAYBACK 

Playback Video Tapes into RHI & PPI 
Separate Sync Pulse From Video 
Playback thru RVP 

Digital Data on Audio, Decode for Analysis 

Display on Front Panel Lights 

Sin/Cos Operate Scopes 

Computer Function to Give 8 or 9 Bit 
Words for PDP-8I 


16 



into 8 or 9-bit words suitable for processing by the PDP-8I computer (an 
integral part of the RVP equipment) for recording on a 9 bit IBM 360 
compatible tape. The sln/cos azimuth and elevation signals with 14 bit 
resolution is outputted on each frame at 150 v for smoothed scope operation. 

The RVP equipment whose functions have been described above 
also had appropriate software written to enable SPANDAR to be operated 
automatically. This permitted the scanning to be made with great economy 
of time. A detailed description of the processor and controller, con- 
structed by Computer Science Corporation, can be found in Appendix B and 
Ref. 21. 

2.4 DISPLAYS 

2.4.1 SPANDAR 

Real-time displays of the radar returns were available at 
SPANDAR. They enabled the investigators to monitor the weather situation, 
to Identify and designate the rainfall areas to be recorded and to serve 
as a guide in later analysis to check cell locations, intensity, area 
covered, etc. Two sets each of Polaroid photographs of the PPI and RHI 
were taken for all sweeps. Examples of PPI and RHI photographs taken 
during the experiment are shown in Figures 2.6 and 2.7. In addition, 
two other scopes were available for the 35 mm camera for backup photo- 
graphs. A photograph of the bank of RHI and PPI scopes is shown in 
Figure 2.8. Log video and sliced (grey scale) video presentations were 

used. Log video showed all rain areas, including the lightest rain. 

Sliced video was used to obtain a quantitative measure of the rain 
intensity and aided in the selection of the more appropriate rain areas 
for studying the individual cells. 

2.4.2 WSR-57 Facsimile Recorder 


The Weather Service of NOAA operates a WSR-57 weather radar at 
the Patuxent River Naval Air Station on the Chesapeake Bay approximately 
90 km N.W. of Wallops Island, Virginia. This radar is part of the 
weather radar network of the Weather Service and is used for weather 
depiction over the Washington, D.C. area. 

The WSR-57 information in the form of a PPI display, is pre- 
sently sent to the Weather Service Station at Wallops Island via a 
leased telephone line and displayed on a WBRR-68 Facsimile Recorder. 
Normally the WBRR-68 transmitter is configured to show a 238 km display. 
This range is quite sufficient to cover the area of interest at Wallops, 
the 140 km radius circle, and to provide adequate forewarning of approach- 
ing weather systems. However, 603 km can be displayed and transmitted 
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Figure 2.6. Example of PPI Photograph During Rain Experiment 
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Figure 2.7. Example of RHI Photograph During Rain Experiment 
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Figure 2.8. Bank of RHIs and PPi? Scopes at SPANDAR 



on request. Figure 2.9 shows the overlapping coverage of the two areas 
of observation. 

One hundred seconds are required to produce a picture on the 
WBRR-68 transmitter. This fixes the highest picture rate at once every 
2 min. A variety of other display rates are available, e.g., once every 
5 min, 10 min, 1 hr, etc. A high quality land line was made available 
from NASA Main Base to SPANDAR. A facsimile recorder was located at 
SPANDAR for use primarily as a forecasting tool. 

A video Integrator and processor (VIP) designed specifically 
to interface with the WSR-57 (not to be confused with the video processor 
designed for SPANDAR) is used at the Patuxent River and other Weather 
Service stations. The VIP is used to provide a contoured display of the 
radar returns which correspond to specified rain intensities. The reso- 
lution element is fairly large, 1.85 km (1 nmi) in range and 2 deg in 
azimuth. A three rpm antenna slew rate is normally used. The range 
displayed when signal contouring is used is from 18.5 km to 230 km. Log 
video is used with the VIP to provide 53 dB of dynamic range. Range 
correction is Included. The contour display is 3 shade, 2 cycle, yield- 
ing six levels of rain intensity up to 10 in. per hr. 

The integrated and contoured data provided a gross check 
against the quantitative data taken at SPANDAR. In addition, the intensi 
ty contour display showed the most intense storms prior to their entering 
the 140 km radius of observation and provided forewarning of approaching 
storm systems. 

3. EXPERIMENT DESCRIPTION 

The first phase of the experiment was carried out during the 
fall of 1972 in which all equipment available at that time was tested. 
Some data were collected during this phase. The full-scale experiment 
was .performed during the summer of 1973 with all of the’ equipment operat- 
ing as described in this report. 

This section discusses the experiments, operational proce- 
dures, calibration techniques and accuracies, theoretical and measured 
(Refs. 10, 11, 13, 24, 25, 26, 27, 28, 29, 30, 31, 32 and 33). 

3. 1 1972 EXPERIMENT 

The preliminary experiment was carried out from 15 September 
to 6 October 1972. Since the Radar Video Processor and Controller were 
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not yet available at that time all storms were recorded on the FR-910 
video tape recorder described in Section 2 and on Polaroid film. 

Frequency diversity components were operating but with a 24-step se- 
quential logic and not with the logic shown in Figure 2.2. 

Data taken during this period were obtained while the radar 
was under manual control. Data obtained in this mode were processed 
"off-line" as described in Section 2.3.1. 

3.2 . 1973 EXPERIMENT 

The 1973 Summer Experiment took place in the period between 9 
May and 31 August 1973. This was the full-scale experiment and the one 
which,; yielded the data analyzed and reported in Volume II. 

< 

From the start of the summer experiment until 29 May recording 
was on video tape. On that date the Radar Video Processor and Controller 
became operational and the experiment went into computer control, or ' 
into the "on-line" processing mode. The processor and control remained 
in operation and has now become an integral part of SPANDAR. Frequency 
diversity for the entire Summer of 1973 was operated in Interlace mode 
as shown in Figure 2.2. 

It was not economically feasible to record every rain within 
the 140' km range of the radar. The philosophy was to collect a smaller 
number but still representative sample of the summer storms of interest. 
The observation times in terms of the number of days, hours of the day 
and frequency of observations were chosen to reflect this approach. 

Two modes of operation were used, as mentioned previously in 
Section 1, each consistent with the two primary objectives. The first 
objective was to obtain a description of the field of cells. Its mode 
was a set of PPI scans ^ usually covering 360 deg, taken at 2 deg eleva- 
tion at one-half hour intervals throughout the time when any rain was 
detectable by the radar within the 140 km range. 

The second objective of this program was to determine the 
detailed spatial and temporal characteristics of single rain cells 
leading to a statistical cell model. In this case, the reflectivity 
structure as a function of altitude and time was measured. This mode 
of operation Involved rapidly probing selected rain cells on a fine 
time scale and with fine spatial resolution. For example, once a cell 
or group of cells was selected the radar probed in stepped elevation 
PPI series through the cell or cells. This is called a PPI raster scan. 
The single-cell mode, then, concentrated fine scale time and space probing 
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in a prescribed, limited volume. This mode was interlaced between the 
statistical mode observations. , 

A full description of the operations plan will be found in 
Reference 13. However, to make this experiment description more mean- 
ingful selected sections of the operations plan which guided the pro- 
cedure in obtaining data are reproduced in Sections 3.2.1 through 3.2.3. 

3.2.1 Equipment Setup 

"A. SPANDAR will be used to collect data from the 1973 simmer 
experiment, “l^e FPS-18 transmitter with frequency diver- 
sity and no paramplifier will be used throughout. Pulse 
length will be 1 psec. Antenna slew rate will be set at 
3 deg sec“^. 

B. The prime data recording will be on digital magnetic 
tape. This equipment, called the RVP, digitizes the 
Incoming video signals, integrates 128 pulses for 
each of the 871 range gates and records the average on 
the DlgiData 1600 bits-per-inch magnetic tape. There 
are two tape transports and data are recorded on each 
sequentially; when one tape is completed the other is 
activated automatically so no data are lost. 

C. The video tape recorder , FR-910, will be in use all 
summer as a backup recording system. At the start of the 
experiment all data will be recorded on the FR-910. This 
requirement may be relaxed after some of the digital 
tapes are processed through the IBM 360 computer and all 
equipment is found to be operating satisfactorily. 

D. An interface system has been designed to allow us to 
operate : 

a. SPANDAR with FR-910 

b. SPANDAR with RVP 

c. FR-910 with RVP 

d. FR-910 with scopes. 

E. Displays and Recording Techniques . Displays of the storm 
cells must be available in real time to facilitate de- 
cisions in the day-to-day operations. The following dis- 
plays will be available: 
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a. Facsimile recorder display of the Patuxent River 
WSR-57 radar operating during the experiment, this 
is called PAXFAX. 

b. There will be two PPI and one RHI scopes with pho- 
tographic recordings. 

c. The video slicer equipment will be Installed at 
SPANDAR to provide quantitative . estimates of storm 
cell intensities. Six slice levels will be used. 
Slice settings will vary depending on the rain 
intensity at the time. Recording will be on Pola- 
roid and 35 mm film." 

3.2.2 Dally Operations: Radar 

"The sequence of events to be followed is (see Ref. 13): 

1. Activate the radar system, scopes, FR-910, RVP inter- 
faces, Fax recorder, slicers, photographic and other 
recording equipment. 

2. Set up RVP and FR-910. 

3. After warm up period calibrate radar, processors and 
recording equipments by stepped attenuation method. 

4. Scan PAXFAX output for indication of rain activity. 

5. Now follow flow diagram shown in Figure 3.1. Tech- 
niques involving scan modes and calibrations on tape 
will be found in Refs. 8, 11, 20 and 21. 

6. After a storm abates the operation will revert to 
the watch mode on PAXFAX and SPANDAR. ^ 

7. Full calibrations will end each day's operations. 

This will include replaying of video tapes to test 
the presence and adequacy of data on tape. Checks 
of Digltapes and RVP operation should also be made. 

8. The interface has controls which must be set before 
each scan. These are required to provide data iden- 
tification of equipment status and experiment param- 
eters. These will include: 
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Figure 3,1. Operation Flow Chart for Rain Statistics Experiment 
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a. 

Year, day, sweep number 


b. 

Time (this is set in automatically) 

c. 

Azimuth and elevation limits 


d. 

Angle change for successive scan 
scan mode 

in a raster 

e. 

Type of run: RHI or PPI 


f. 

Number of scans in a sequence 


g- 

Radar status: prf, pulse width, 

data validity 

scan rate. 

h. 

Calibration levels 



9. A detailed written operations log will be kept to record 
the sequence of events during a dally experiment. These 
logs will Include all the information a. through h. above 
as well as supplementary Information as equipment mal- 
function, if any. Standard RHI and PPI logs will be kept 
for Polaroid and 35 mm photographs. The FR-910 also has 
an audio channel on which data a. through h. -above should 
be voice-recorded. 

10. Data handling will require care since there will be so 
much of it and in such a variety. Video tapes will 
always be documented and kept at the SPANDAR site since 
it requires the FR-910 for playback. Dlgitapes will be 
transferred SAP to APL for processing. Polaroid photo- 
graphs, logs, weather data, strip chart records calibra- 
tions and weather notes will be documented and returned 
to the RAP project at APL once per week. 

11. • The SPANDAR operation, starting May 7, 1973, will be on a 

3:00 p.m. to 11:30 p.m. shift. Four operators will be 
assigned to that shift on a rotation basis. The Site 
Manager (SM) will be on a shift which is split between the 
rain and the clear air research. The Project Director 
(PD) will start his daily shift with the 2:00 p.m. weather 
briefing." 


27 



3.2.3 


Graphical Aids 


"A. Figure 3.2 is a height-range chart for various antenna 
elevation angles. Range is in nmi out to 100 and height 
extends up to 70,000 ft. The chart was drawn for a 4/3 
earth radius. 

B. Radar power vs reflectivity factor. This chart is Figure 
3.3. On the upper chart we have radar power received, 
dBm plotted against reflectivity factor, dBZ for selected 
values of distance from 10 to 140 km. Note these are for 
nominal values of x and k. On the lower chart we have 
rain rate vs reflectivity factor for selected Z vs R 
relationships. Ordinates are given in mm hr“^ and in 
hr" ^ . 

C. Figure 3.4 is included to aid in normalizing power values 
for the various ranges. Power received in dBm at 1 nmi 
is plotted against range in nmi for various values of 
radar received power over the range of Interest for this 
summer’s experiment. A -30 dBm signal at 10 nmi for 
example, is found on this graph to be -50 dBm (normalized 
to 1 mi); this corresponds to a -13 dBm signal at 70 mi. 

D. Figure 3.5 is plotted to give a quick look at the change 
in elevation angle required to give an approximate change 
in vertical height of 1 km and 2 km. These are approxi- 
mate because the angle change for a given height change 
is different at different altitudes at the same slant 
range. " 

Note that the flow diagram. Figure 3.1, contains a Z-R mode of operation 
Although efforts were made to use a disdrometer and a rain gage along 
with the radar to obtain a Z-R relationship for various rainstorm in- 
tensities, the disdrometer under loan from another organization, did not 
operate satisfactorily. That portion of the flow diagram which is 
shaded was not used in any significant way; therefore it may be con- 
sidered "inoperative" for this experiment. 

3.2.4 Meteorological Instrumentation 

Two different types of meteorological data were required in 
support of this experimental program. One set of data concerned the 
mesoscale and synoptic features of the weather system and the other set 
of data the point measurement of the rainfall. In the former case, the 
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Figure 3.2. Height Range Chart for Various Elevation Angles 
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Figure 3.3. Radar Power Versus Reflectivity Factor for Various Ranges (upper chart); 
Rain-Rate Versus Reflectivity Factor for Various Z-R Relationships 
( lower chart) 
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Figure 3.5. Elevation Angle Change to Produce A Change in 
Altitude As A Function of Range 
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facilities for collecting and recording the data were in existence and 
available at Wallops Island and other Weather Service stations in the 
area. In the latter case, the instrumentation was installed and opera- 
ted by APL personnel. Since the two sets of data serve entirely differ- 
ent purposes they will be discussed separately. 

Mesoscale and synoptic weather data were required for three 
reasons: (1) for the purpose of forecasting; (2) for the purpose of 

relating the radar data and the point rain measurements to the broad 
meteorological development; and (3) for use in categorizing the rain 
types and in extrapolating the results to other regions. 

The Weather Service station at Wallops was used as the primary 
source of weather information. Surface observations, radiosonde data 
and surface measurements were made available along with all the facsimile 
charts from the Weather Service network. In addition, weather briefings 
were held each morning which included local forecasts. 

Surface observations and radiosonde data from other Weather 
Service stations such as Washington, Hatteras and Greensboro were also 
assembled. The radar data from the Weather Service radar at' Patuxent 
River was available at the Wallops Island Weather Service station and 
at SPANDAR. 

3.2.5 Point Measurements 


The radar-rain relationship, commonly referred to as the Z-R 
relation, is needed to translate the radar reflectivity (Z) measurements 
within the storm cells into rain rate (R) and total water content (M) . 
Rain rate alone, as determined from ground-based rain gages such as 
tipping buckets, is not sufficient. Joss et. al . Ref. 34, has found 
good agreement when the drop size distribution, rain rate and radar 
reflectivity were categorized in terms of the type or rain, l.e., driz- 
zle, widespread and thundershower. Categorization is needed if extrap- 
olation to other areas is to be attempted. 

Point rainfall measurements were made and were an Integral 
part of the operation. A tipping bucket rain gage was installed at the 
VOR site, about 12 mi north of the radar site. These data are augmented 
by rainfall data from other existing weather stations, l.e.. Wallops 
Station, Patuxent, Dulles, Atlantic' City and Richmond. As part of the 
experiment plan specific instructions were issued to the Project Director 
related to the meteorological requirements. Some of these instruction^ 
are quoted below* 
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3.2.6 Dally Operations ; Meteorological 

"A typical sequence of events Involving the PD and meteoro- 
logical aspects of this experiment Is as follows: 

"1. The PD visits the forecaster on duty at about 1400, 

before start of the 1500 shift. The PD obtains Infor- 
mation on whether rain/showers are expected within 140 
km of SPANDAR for the period 1500 to 2330, and. If so, 
from what direction It Is expected to first appear. If 
forecaster states the type of rain expected, this should 
be recorded on the metlog. Also, the forecaster should 
be asked If he expects clrrlform clouds' (but not anvil 
cirrus), and his answer Is recorded. If It Is already 
raining at 1415, or has been raining, find the likelihood 
of another rain period before 2300 EST. This Information 
Is desired for guidance to the PD and need not be recorded. 

2. The PD will use the Information transmitted by the Patuxent 

radar to our site. The transmitted copies of PAXFAX 
should be retained for archiving. If date or time Is not 
clearly legible, the PD should enter these on the copy. 
Note: The Patuxent radar information Is advisory. We 

may be In better position than Patuxent to observe rain 
first appearing from the South to East. There may be 
convective situations In which showers develop spontan- 
eously In our observation region so there Is no early 
warning from the Patuxent radar. 

3. When we are actually observing rains or showers with 
SPANDAR, the PD informs the Weather Service that we are 
observing precipitation. This is needed because this 
initiates a special set of information from them, re- 
quired for our analyses. 

4. The PD (possibly requiring assistance of site person- 
nel) has the responsibility for maintenance of the rain 
gage. He visits the gage to change charts, check tipping 
bucket operation, general operation of the gage measure- 
ment and other checks described below." 

"a. Chart change. Prior to removal of old chart he jogs 
pen, to show record termihatioh, and writes date and 
time (EST) on the chart. Correspondingly on new 
chart, pen jog denotes record start; date and time 
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(EST) are entered on chart. Chart is moved to cor- ; 
respond to the correct time, allowing for backlash. 
The box containing the chart record is marked with 
location, "rain gage," date and time of beginning , , ' 
and end of record and rain amount (see below) . 

While a chart will easily record for one week it is 
desired that the chart be changed during a non-rain 
period, or if unavoidable (say, near end of chart) 
during a light rain since we want to avoid loss of 
any rain tip information. It is desired to obtain 
rain gage records through the period of the summer 
experiment, whether or not there are radar observa- 
tions. 

b. Miscellaneous. The inking of the charts are verified 

from examining the end of the old chart . Fresh inks 
are added to pens as necessary. Pens are replaced 
and ink flow checked as necessary. Determine, if 
possible, if there has been a power interruption, 
and, if so, if the time of power outage is known or 
can be estimated. A note containing this informa- 
tion should be included in the met log. The tipping 
bucket is tipped manually 10 times at the end of the 
chart record. If the chart does not show 10 tips, 
attempt to repair or report problem. If there is 
water in the tipping bucket, estimate if it is less 
than half full in which case it is neglected or more 

than half full in which case note at the end of the 

chart that there is a residual .01 in. rain. How- 
ever, the manual tipping of the bucket is not done 

if it is raining. 

c. Rain Amount, The total rain that accumulated during 

the period of chart recording — or more often — is 
measured. The use of a standard Weather Service 

inner cylinder of the 8 in. gage, or the use of a 

graduated glass cylinder, is recommended. This 
information, either as a total for the entire chart 
or as several separate readings are entered at the 
end of the applicable chart and on the chart box. 
After measurement, the water is dumped. 

d. ■ Site Examination. The gage should stand vertical 

with opening horizontal. The windshield around the 
gage should show no sag. If the site owner wants 
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the grass mowed in gage area. It is preferred that 
we do it to be sure no damage is done to the equipment. 

e. Visits and work done are recorded on the metlog." 

"5. The PD will attend a special 1400 weather briefing at the 
Wallops main station. He will be advised of the likeli- 
hood of showers during the afternoon and evening of that 
day. He will study the latest synoptic chart, radiosonde 
and other aids to help guide the experiment for that day. 

6. He will proceed to SPANDAR to organize the operation. 

His decisions from this point onward will be determined 
by continuing study of the PAXFAX, SPAHDAR PPI Polaroid: 
sliced photos showing the progression of the storms and 
the sliced video Polarolds along with continuous contact 
with the Wallops station weather forecaster. SPANDAR PPI 
scans at a 2 deg elevation angle will be made at least 
every half hour when precipitation is predicted. 

7. The PD receives various materials from the National 
Weather Service Office, Wallops Island, as indicated in 
the Facilities request. If the material is hardly legible, 
or missing, he should attempt to get suitable copy from 
the Weather Service. Checks should be made that each 
Item is dated. Each day's records should be put in 
separate folders. 

8. The PD will maintain the metlog. A form will be provided 
for this purpose. While all indicated items on the form 
would be complete, entries on where and when rain or 
showers are observed on the radar should also be Included 
plus whatever the PD regards as significant, of special 
Interest, etc. The log plus other material of each day 
should be placed in dated folders and checked that all 
the material has been collected into the folder. Chart 
rolls, film and cassettes are kept separately and with 
the folders are brought to APL by the PD." 

Other aspects of the meteorology program will be discussed in 



3.3 


RESOLUTION, ACCURACY AND CALIBRATION 


This section describes features of the equipment relating to 
resolution and accuracy. Considerable emphasis was placed on these 
aspects since the goal was to have an overall system accuracy of better 
than + 5 dB. At the end of this section it is shown that calibration, 
stability and processing errors were + 4.05 dB peak and less than 1.65 
dB standard deviation for this experiment. 

3.3.1 Spatial Resolution 

The overriding consideration in the selection of the SPANDAR 
radar for this program and the processing and recording techniques was 
the preservation and exploitation of the high spatial resolution avail- 
able. The spatial resolution cell is defined to be one radar beamwldth 
in azimuth or elevation angle by one pulse length in range. The radar 
calibration included detailed beam pattern measurements so that the 
actual beamwldth is known. Ref. 32. The nominal width of the radar beam 
along with the aspect ratio of the spatial resolution cell as a function 
of range are shown in Figure 3.6. The 1 ps pulse length was chosen to 
take full advantage of the beam dimensions at close range. Note that 
the minimum observable range was 10 km due to ground clutter. At long 
ranges the beam becomes quite wide, e.g., around 1 km at 140 km range 
and the aspect ratio is also large. The effect of this beam-broadening 
on the accuracy of the reflectivity measurement is discussed later. 

3.3.2 Dynamic Range 

In this program it was critical to have the capability of 
measuring the very large reflectivities associated with hail-producing 
thunderstorms. The largest reflectivity values recorded are near 10^ 
mm®/m^ (Donaldson and Tear, Ref. 35). The radar and processor had to 
be able to measure Z values of this magnitude at close range to ensure 
agklnst the loss of a rare but extremely important event. On the other 
hand, the dynamic range of the radar also had to be wide enough to allow 
detection of less Intense rain at 140 km. 

The system dynamic range was determined by the log IF ampli- 
fiers. As mentioned earlier, the extended range log IF amplifier had a 
linear range of 90 dBm. 

/ 

The variation in power received as a function of range and the 
equivalent reflectivity factor was shown in Figure 3.3. The radar ' 
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constant, k, and the antenna beamwldth, 0, are nominal values. The 
actual values were determined by the calibration program described 
later. Also included In Figure 3.3 is a plot of equivalent reflectivity 
factor vs rain rate for several measured Z-R relationships. The figure 
shows the range of the different parameters which can be measured in the 
manner described below. 

A conservative estimate of the radar noise level is about -98 
dBm (without the paramp) . This assumes a noise figure of 10 and a 4 MHz 
bandwidth throughout the system. The true noise figure was determined 
by calibration during the experiment. Setting the lower limit of the 
linear range of the IF amplifier at a power received of -95 dBm just 
above the system noise, and allowing 10 dBm for fluctuations in the 
received signal, the lowest average power received which could be measured 
with confidence is -85 dBm. At the 140 km range, then, the minimum 
reflectivity factor is 27.5 dB. If one assumes, for example, the Joss 
model for thundershowers, the minimum rain rate which could be measured 
at the maximum range is approximately 1 mm hr"^. 

The upper, limit of the 90 dBm linear range of the IF amplifier 
was at -5 dBm. Again allowing 10 dBm for signal fluctuation, the maximum 
average power received that could be measured with confidence was -15 
dBm. This corresponds to a maximum measurable reflectivity of 75 dB at 
the closest range for which data were recorded, 5.4 nml (10 km). Thus, 
the 90 dBm dynamic range of the log IF amplifiers covered all the reflec- 
tivity factors and rain rates of interest. 

3.3.3 Measurement Resolution 

The incoming radar video is first digitized to 6 bits by the 

A/D ^converter in the on-line data processing system. For a 90 dBm range 

of signal levels (discussed in the previous section) , the quantization 

step size is 1.4 dBm (=90 dBm/64). The output of the integrator is an 8 

bit word so that the output resolution is 0.35 dBm (=90 dBm/256). The 
uncertainty in the integrated power received due to the above trun- 
cations or quantizations is discussed in the section that follows. 

3.3.4 Data Accuracy and Precision 

The various sources of error contributing to the overall ac- 
curacy and precision of the radar measurements are discussed in the 
sections that follow. The estimated individual errors are worst case 
estimates and are not likely to be exceeded. Their sum is, therefore, a 
most conservative error estimate since individual statistical errors 
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will tend to cancel to some extent and systematic errors may be accounted 
for in data reduction. 

3.3.5 Statistical Error 


An extensive discussion on the probability distributions of 
fluctuating echoes from randomly distributed scatters has been given by 
Marshall and Hitschfeld (Ref. 16). Only the pertinent results will be 
given here. 


Fluctuations occur in meteorological targets having constant 
reflectivity because the reflectors move with respect to each other and 
the phases of their backseat tered signals change in a random fashion. 

An estimate of reflectivity can be obtained by averaging a number of 
independent echoes. The accuracy of estimated echo power, Jj^^, depends 
on the number of samples included in the average. The distribution of 
the average k samples, around the true value, is given in Figure 
3.7 taken from Ref. 16. These curves were derived analytically for A^. 
However using a Monte Carl o tech nique they show that these same distribu- 
tions are appropriate for log A^ which is too complicated to handle 
analytically. Figure 3.8, derived from Figure 3.7 gives 50%, 80% and 
95% confidence limits for the average k samples of log A^ to be within e 
of the true value. For example, an average of 100 samples of log A^ 
falls within + .4 dB of log A^, minus a known constant, 80% of the time. 
Note that the accuracy improves very slowly for values of k greater than 
roughly 60. The RVP system integrates 128 independent pulses in the 
average received power that is recorded. Thus, the statistical error in 
the recorded data from this source should be slightly less than the + .4 
dB quoted above. 


3.3.6 Quantization Error 

A digital video integrator was used in the RVP system. There- 
fore, the errors introduced by quantizat ion of both the integrator in- 
put, log A^, and the integrator output, log A^, must be estimated. 

Austin and Shaffner (Ref. 36) have considered the effect of 
quantizing the input to a digital integrator. It is shown that the 
logarithmically quantized average depends on the width of the quantizing 
intervals, and to a much lesser extent, the position of the true mean 
intensity within an Interval. Thus they arrive at a relationship between 
input quantization error and the width of the intensity level interval. 
Figure 3.9 shows this relationship. Note that the error for no quantiza- 
tion whatsoever is -2.5 dB which agrees with Marsh all and Hitschfeld *s 
result (Ref. 16) for the difference between log A^ and log A^ for analog 
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Figure 3.9, Error Due to Input Quantization 

averaging. The assumption implicit in this discussion is that all 
statistical fluctuations have been smoothed by including a very large 
number of samples in the average. Therefore, if the appropriate correc-^ 
tlon is selected from Figure 3.9, the error from the quantization of log': 
at the input to the integrator can be corrected completely. For the ■ 
1.4 dBm quantization step size in the A/D converters the input quantiza- 
tion error (or correction) is 3.2 dBm. > : 

The maximum error contribution from the quantization of the 
integrator output is simply the value of the least significant bit. The 
output of the integrator was quantized into 256 levels over a dynamic '' 
range of 90 dBm so that the output quantization and, thus, the maximum' 
error was 0.35 dBm. 

3.3.7 Error Due to Reflectivity Gradient 

An additional error is introduced when samples from a log 
receiver are averaged when the radar sample volume contains* a reflec- 
tivity gradient. The conventional probability distributions for scat- 
terers uniformly spread across the pulse volume (Ref. 16) are distorted 
when such gradients are present. Both the variance and the mean are ' 
affected in the case of log averaging. 
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This effect has been analyzed by Rogers (Ref. 37) for two 
different reflectivity gradient models; model 1 assumes a linear varia- 
tion In log I across the pulse volvane and model 2 assumes a linear 
variation in I. Model 1 is believed to be the most realistic. The 
relevant results of this analysis are presented In Figure 3 . 10 which 
shows how the estimates of variance and mean depend on the gradient 
parameter k. This parameter Is the ratio of the power received from the 
most strongly reflecting region in the pulse volume, I 2 , to the power 
from the weakest reflecting portion, Ij. Other parameters in the figure 
are defined as follows: 


J 


<J> 



the average of log I for scatterers having 
constant reflectivity (equal to I^, an average 
across the gradient) . 

1 

the' estimate of log I obtained from an average of 
log I when a gradient Is present across the pulse 
volume. 

signal variance for scatterers having constant 
ref lectlvlty <I^> . 

the measured signal variance of log I when a 
reflectivity gradient is present. 


Gradients of 20 dB/km are not uncommon in thunderstorms. Such 
a gradient located at 140 km from SPAMDAR represents a vrorst case situa- 
tion. The beamwldth at that range is about 1 km so that a 2D_ dB varia- 
tion in reflectivity will be used in this example. -T-he-cuirves evaluated 
at the corresponding k value, 100, shows that a bias of 3.-2-dB_is intro- 
duced, and the variance is Increased by about_3 dB for the more repre- 
sentative model 1. Thus in this worst case one would underestimate the 
Z value by about 3 dB, and the statistical uncertainty of~~. 4 dB for the 
uniform distribution of reflectors Increases to about .8 dB or less for 
80% of the time. Note that these errors improve significantly at dis- 
tances of about 65 km where k reduces to 10 for a 20 dB nml~^ gradient. 


Although the foregoing shows that errors of 3 dB or less are 
possible in regions of high gradient, the Z value of most interest, 
l.e., its maximum value, will be least affected by this type of error. 

3.3.8 Integration Error 

I An error contribution is produced by the integration of N 

Independent pulses as the beam sweeps through the target. If we assume 




a true averaging (rectangular weighting) digital integrator, the output 
of the integrator is 


N-1 




where is the i^^ input to the integrator in dBZ. Donaldson's (Ref. 
56) storm profile for model 1 is of the form 


e = dBZ„j-mX2 

where dBZ„ is the dBZ value in the storm core, m is the constant de- 
m 

fining the apparent storm size as observed with a finite beam, and X is 
the distance from the storm center. The i^^ input to the integrator 
will be the sum of the mean Z value at a distance iojRT from the thunder- 
storm core plus a noise component, n . With no loss of generality the 
assumption is made that the beam is directed at the storm center for the 
sample 1 = o. In the above, oi is the scan rate in rad/sec, R is the 
range, and T is the time between independent samples. Thus, the i*^^ 
input to the digital Integrator becomes 

« dBZjjj - m(X-iu)RT)2 + n^^ 


so that 


since 


_ N-1 

^0 = " I 2 

i=o 


, N-1 

N 2 "i 
i=o 


0 . 


Carrying out the sum produces the result. 


S = ‘^^^m - l" 2 


[X - i coRT(N-l)l^_l m(o)RT)^ (N^-l) 
12 J 12 
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which shows a bias error of 


= - ^ m(o)RT)2 (N^-l) 

and a displacement error of 

6X = 1 ojRT(N-l) . 

2 

Worst case values for the present experiment are as follows: 

■ ■ 0 ) = 360°/92 sec = .07 rad/sec 

R = 140 km 

T = 10“^ sec 

N = 128 

m % 1 dB/km^ (Ref. 56). 

These values lead to 

6e ~.2 dB, 

0 

and ■ ' 

6X ~ . 6 km, 

both of which were acceptable for the present experiment. 

3.3.9 Beam Shape Error 

Extensive work has been done by Donaldson and Tear (Ref. 35) 
and- Donaldson (Ref. 56) in estimating the effect of smoothing by the 
radar beam. Functions representing seven storm profiles were convolved 
with an idealized beam pattern to provide the smoothed profiles used to 
estimate distortions in cloud profiles and errors in Z measurements. 

Extrapolations from the results of Donaldson and Tear (Ref. 

35) indicate that a .5° beam would contribute errors in the measurements 
of maximum Z value of less than 2 dB at 100 km range and less than 4 dB 
at 200 km. The seven models used were typically about 4 or 5 km in size 
at their 3 dB points. Therefore, it is believed that variations in cell 
reflectivity across the beam only infrequently result in Z^^^ errors 
greater than 2 to 3 dB in the vicinity of 150 km. 
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Severe distortions can result in the measured storm profiles 
at distant ranges because of the finite beamwldth. Results for Donald- 
son's model 1 (Ref. 55), one of his most realistic, are shown in Figure 
3.11. The figure shows the effect at longer ranges of the side lobes 
illuminating the intense core when the main beam is well above the storm 
top. For this reason, accurate measurements of storm height should be 
performed only at a ;close range and with a narrow beam radar. Note that 
for this model contours as much as 5 dB below the measured value are 

not strongly dependent on range. 

The beamwidth for SPANDAR is 1.0 km across the 3 dB points at 
140 km range. This is the maximum acceptable range that would enable a 
rather small cell of about 3.5 km diameter at points 3 dB down from to 
be measured with acceptable resolution and accuracy. 

3.3.10 Radar Stability and Accuracy 

The nature of the rain cell statistics program in which radar 
reflectivity is measured and transformed into rain intensity levels 
required that every possible effort be expended in calibrating the, radar 
sensor for absolute measurements. This requirement for accuracy of 
measurement is commensurate with the fine spatial resolution available 
at SPANDAR with its 0.4 deg beam and 1 ys pulse length. 

Since measurements were taken over an extended period of'^time, 
the drift and stability of the sensor (the radar system) was important 
both on a short-term (hours) and a long-term (months) basis. Observa- 
tions taken during a given sequence of scans had to be accurately relat- 
able to each other and to other series taken several months apart. 

The calibration program for the SPANDAR radar is outlined in 
Table 3.1. There were two major elements in the program, the pre- 
experiment program and the during-experiment program. The pre-experi- 
ment program was primarily concerned with the measurement of system 
cl^kracteristics, e.g., antenna patterns, gain, transmitter bandwidth, 
p0wer, transmission losses, receiver noise figure, bandwidth, etc. The 
details of this program are contained in Vann (Ref. 32). Preliminary 
results are in Vann (Ref. 27). 

In establishing an accuracy figure for the radar, attention 
was given to the following: 

a. Location at which the calibrating signal is 
injected ^ 

,b. Location where the measured gain is referred 
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ANTENNA SMOOTHED CORE REFLECTIVITY (dB) 



0 2 4 6 8 10 12 


DISTANCE FROM STORM CENTER (X,nmi) __ 

Figure 3.1 1. Effect of Beam Smoothing on Donaldson's CloujlModel 1 (Ref. 56) 
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Table 3.1 


r 

Calibration 


A. 

i 

INITIAL MEASUREMENTS > 


1. 

Confirm Linearity of Polarization 1 


2. 

Measure Line Losses Including Rotary i Joints 


3. 

Measure Antenna Gain and Pattern ' 


4. 

Measure Receiver Gain and Bandwidth i 


5. 

Design and Fabricate a Method for Finding 



Scanner Position When Stopped ' 

B. 

PERIODIC 

CALIBRATIONS ‘ 


1. 

Use Boresight Tower 1 


2. 

6 in. Sphere 1 

C. 

DAILY CALIBRATIONS AND TESTS 


1. 

Measure Minimum Detectible Signal j 


2. 

Monitor Transmitted Power and Pulse Shape 


3. 

Set Up Log Strip j 


4. 

FR-910 Check for Optimum Operation ! 


5. 

Test RVP; CW Signal with Ramp to Test All Bins 



and All Steps in Each Bin; Use Autbmatic 
Stepped Attenuator 

D. 

CALIBRATION AND TESTS FOR EACH RUN ! 


1. 

Transmitted Power ! 


2. 

Frequency Diversity Check: Deviation and Power 

E. 

CONTINUOUS CALIBRATIONS AND CHECKS 


1. 

Transmitted Power 


2. 

High Level Signal ' 


3. 

Low Level Signal ' 


4. 

Receiver Noise (Scope Monitor) 

i 


! 
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c . Accuracy of the gain measurement 

d. Accuracy of the measurement of insertion loss • , 

between the locations where the gain is measured 

and ithe calibrating signal is injected 

e. Stability of the transmitted power and the 
accuracy of the calibrating signal. 

The calibrating signal was a sample from the transmitter input., 
and the accuracy of the losses in the couplers and associated lines had 
to be measured. Since this calibrating signal was injected into the 
mixer input of the receiver, gain instabilities and losses in the receiv- 
ers Influenced the incoming signals in the same way as the calibrating • 
signal. Hence, the characterization of the receiver properties played a 
minor role in the accuracy question. 

' • 1 ' • 

Since frequency diversity was used, the variation in trans- ■ 
mitted power as a function of frequency had to be measured. The VA-57 
klystrons used in the FPS-18 transmitters could be tuned to be uniform 
in power to within 3 dB over a 25 MHz frequency range. The slight 
reduction in average transmitted power due to frequency diversity was., 
taken into account in system calibration. 

The second part of the calibration program concerned the moni- 
toring and/or logging of critical points in the radar system during the , 
experiment. Certain parameters such as transmitted power, calibration , 
levels, and receiver noise were monitored and recorded continuously. -,i , 
Other items were checked and/or logged less frequently such as, daily or ■ 
weekly. The details of these measurements are contained in Ref. 32.,. 

3.3.11 Calibration 


The calibration program comprised a number of measurements 
which may be divided into the following major phases: 

Phase 1: Pre-program measurements performed 

before the rain measurement program 
was initiated. 

Phase 2: Measurements performed before, during, , .. 

and after each rain measurement day. 

In the following paragraphs, a description of the types of measurements, • 
and the procedures followed are given for each of the phases and ref- . , ; 
erences are given to applicable memoranda characterizing the associated 
measurements . 
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3.3.11.1 Phase 1 - Pre-Program Measurements . The aims of this phase 
were to establish the operating characteristics of SPANDAR as well as to 
enable a subsequent system of monitoring to achieve a high level of 
accuracy during the experimental period. A complete r’eport on these 
efforts is given in Ref. 32 and is included in this Final Report as 
Appendix C, Volume III. , 

Specifically, a large number of measurements were made cover- 
ing the transmitter, antenna, and receiver along with the associated 
waveguide transmission lines, monitoring couplers, coax lines and at- 
tenuators. The measurements were performed initially to determine the 
general condition of the microwave system and to identify any component 
that had to be replaced due to high VSWR or insertion loss. The meas- 
urements were then repeated after replacing several sections of wave- 
guide. 


Some concern had been expressed regarding the second and third 
time around rejection of radar returns due to the frequency separation 
between adjacent pulses. In order to provide this Information the radar 
preamp response was measured and the entire receiver chain was measured 
over an extended frequency range at signal levels covering the full 
dynamic range of the system. A subsequent analysis (Ref. 38) demonstrat- 
ed that second-time-around errors could be neglected. 

The FD chassis was a key element in the experiment and a large 
number of tests were performed on the unit in order to determine th'e 
stability, frequency linearity, power output variations across the jfre- 
quency band and the stability of the timing system (Ref. 14). Some 
problems were observed in the timing system stability and the circuitry 
was modified to a counter-controlled delay generator with a much improved 
stability. 

I 

Extensive antenna patterns were measured covering the principal 
planes at several frequencies, along with backlobe and spillover lobe 
response and a set of stepped elevation patterns were obtained from 
which antenna contours were generated. Antenna gain measurements have 
been reported in Ref. 30. 

A system of continuous monitoring and calibration was devised 
for the SPANDAR radar which consisted of continuous measurement and 
recording of certain par£imeters, periodic measurement of other parameters, 
periodic injection of a known level signal from a source external to the 
radar, continuous injection of noise samples and samples of the transmit 
frequency timed to occur in the radar dead times. These are alluded to 
in the description of Phase 2. 
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3.3.11.2 Phase 2 - Measurement Program During Experiment . The calibra- 
tion, program followed is essentially that outlined in the operations 
plan described in the Operations Test Plan, Ref. 13. The FD unit rep- 
resented a vital component of the system and hence the continuous moni- 
toring of this unit was considered essential for guaranteeing the validi- 
ty of the experimental results. 

Figure 3.12 is a block diagram, which among its other features 
describes the various outputs of the FD chassis and how this unit re- 
lates to the other components of the radar system. Also shown in this 
diagram are the nominal power levels and pulse widths at the various 
outputs of the FD chassis and at the other stages of the system. The 
major procedures and the associated calibration results are listed in 
the following paragraphs where references are made to the numbered 
components [ ] of iPigure 3.12. 

I 

1. Before each rain measurement day, a known test signal was 
injected daily into a standard gain horn on the 250 ft met- 
eorological tower at Wallops Island and the received power at 
the mixer input was monitored. Such measurement represented a 
daily check of possible fluctuations of the antenna gain and 
line losses down to the mixer input. The fluctuations observed 
during/ the entire summer varied peak to peak within +0.3 dB 

of a given average level. A description of this test is given 
in Ref . 30 . 

2. Step calibrations of the receiver-processor system were perform- 
ed in 3 dB steps over a 90 dB dynamic range (-100 to -10 dBm) 
before, during, and after each rain measurement day. The 
outputs, in digital form, were monitored in real time on the 
teletype and simultaneously plotted to establish the existence' 
of possible anomalous effects or deviations from logarithmic 
linearity. On occasion, when such anomalies were observed, 
they were subsequently corrected before continuing with the 
experiment. The system linearity was found to be well within 

+ 1 dB throughout the entire season. 

3. The forward and reflected power, [s ] and [&], as well as the 
High Cal. Monitor, [9], were continuously monitored with a 
strip chart recorder during the entire rain measurement day. 

The transmitted power was found to fluctuate within +0.1 dB 
relative to a peak power of 1 MW (as referenced at the power; 
output monitor). The calibrate output fluctuated daily within 
+ 0.2 dB, although a long-term drift was observed. This drift 
was accounted for by noting its daily level and adjusting the 
FD unit to give a constant calibrate output. 
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Figure 3 2/ Block Oiagram of Pertinent Gomponents of SPANDAR with FD Chassis 
















4. 


During each one-half hour interval the following levels were 


logged. 



a. 

Local oscillator crystal currents 



b. 

Forward and reflected power [5] and 

b] 


c. 

FD chassis LO monitor output [s] 



d. 

FD chassis High Calibrate Monitor output 

[9] 

e. 

FD Transmit Monitor [?] 



f. 

High and Low calibrate video outputs 

(output of [3] ) 

g- 

Center frequency, f^. 




i 


The above records were; ■ examined daily for diagnostic purposes to estab- 
lish the extent of any deviations from the nominal levels. In general, 
all levels remained sufficiently close to their nominal values during 
the entire rain program. 

' I 

3.3.11.3 Calibrate and Noise Signals Recorded ori Tape . Averages of 128 
samples of the high and low calibrate signal as well as receiver noise 
were recorded continuously on digital tape in defined bins associated 
with each transmitted pulse. The high calibrate signal is at a frequency, 
f , with the leading edge occurring at a time, t^ + 1002 ps (t^ = time 

or leading edge of transmitted pulse) , the pulse width is 4 ps and the 
pulse power is 0 dBm (assuming zero attenuation at [z ] . The low cali- 
brate signal is a noise pulse whose leading edge occurs at a time, t^ + 
1020 ps, has a pulse width of 15 ps, and a power level of -70 dBm (aver- 
age power over a pulse). The output is at the same location as that of 
the high calibrate pulse [ lo]. The combined signals are injected into 
the receiver mixer through a 20 dB coupler [ll] resulting in power 
levels of -20 dBm and -90 dBm, respectively, at the mixer input terminals. 

These averaged high and low calibrate signals were used in 
subsequent programs to establish intermediate power levels through in- 
terpolation procedures. This calibration scheme allowed for the com- 
pensation of either short- or long-term receiver gain changes. 

3.3.11.4 Radar Calibration Constant and SPANDAR Error Budget . In Ref. 

30 are described the measured parameters of SPANDAR which appear in the 
radar constant. In addition, the total error budget is derived and the 
uncertainty in the reflectivity factor, Z, established. Although these 



error uncertainities were in part based on initial measurements during 
the early stages of the rain statistics program, subsequent measurements 
resulted in uncertainties consistent with the earlier levels. 

The resultant reflectivity factor as extracted from the radar 
equation was shown to follow the relation, ' ^ 

Z(dB) = (dBm) + 20 Log^Q r(nmi) + (1) 


with. 


Ki = 76.3 dB 


( 2 ) 


where (dBm) represents ten times the logarithm of the average power 
expressed in dBm. This quantity was extracted from the average 'of the 
logarithmic power through the addition of a 2.5 dB correction term. 

An additional 0.9 dB has also been added to correct for; quan- 
tization errors, see Ref. 30. in (1) represents the resultant! radar 
constant which accounts for the correction in Pj-j^, transmitted power, 
gain, receive and transmit losses, wavelength, beamwidths, rain refrac- 
tive index, pulse width and any associated constants, r(nmi) represents 
the range (expressed in nautical miles) and represents the two-way 
atmospheric absorption primarily due to collision broadened resonances 
of the oxygen molecules. This quantity is dependent upon the range of 
the target, and elevation angle of the beam. This term may give levels 
as high as 2 dB (for grazing angles at 140 km range) at S-band and was 
accounted for using a CRPL exponential reference atmospheric model 
(Ref. 28). 

3.3.12 Error Summary 

In Table 3.2 is given a summary of the error budget (extracted 
from Ref. 30). 


The first column lists the peak error levels and the second 
the upper limit values for the standard deviations. The radar calibra- 
tion stability errors Include the uncertainties in the gain, line losses, 
beamwidths and transmitted power. The processing errors include the un- 
certainties in the system linearity, output quantization, statistical 
averaging and interpolation due to instabilities of the high- and low- 
calibrate signals. The phenomenological errors include: beam smoothing 

(associated with the smoothing by the nonuniform beam pattern of the 
reflectivity profile); Integration errors (associated with the inte- 
gration of 128 pulse occurring simultaneously as the beam sweeps through 
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nonuniform scattering regions) ; reflectivity gradient errors (associated 
with the averaging of the logarithm of the power for cases in which the 
radar pulse volume contains a reflectivity gradient) ; and the absorption 
error (associated with the uncertainty on the modeling of the O 2 correc- 
tion term) . ^ 


Table 3.2 
Error Summary 



Peak 

S.D. 

1. Radar calibration and stability 
errors 

2.0 

1.0 

2 . Processing errors 

1 

2.05 

1.3 

Total peak of (1) and (2) 

4.05 


Resultant Standard Deviation 
(1) and (2) 


1.65 

3. Phenomenological and other 
uncertainties 

-2.5 + 0.2 



4. METEOROLOGICAL CONSIDERATIONS 

I 

Weather entered the rain statistics program in several ways. 

It was of prime consideration in carrying out the experiment properly, 
i.e., the Project Director needed accurate forecasts for planning and he 
needed a knowledge of the present weather pattern during daily opera- 
tions. Then, of course, he required measurements of the presence of, and 
amounts of, rain during the course of radar data accumulation. In 
addition, for analytical purposes after the experiment, the analyst 
needed temperature, moisture, wind and precipitation information in 
three dimensions and changes of these parameters with time. With these 
the analyst could categorize the rain structure to carry out the impor- 
tant task of extending the statistical measures of rain cells to other 
climatological regions . 

This section discusses some relevant features of precipita- 
tion, their characterization and what was needed to carry out this 
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effort in the present experiment. More detail along those lines will be 
found in References 39, 40, 41 and 42. 

4.1 RAIN CHARACTERIZATION 

There were many rain characterization schemes in use or sug- 
gested. For example, Austin and Houze (Ref. 43) used the division: 
large mesoscale areas, small mesoscale areas and cell. Austin (Ref. 44) 
used thunderstorm rain and nonthunderstorm rain for characterization 
(and includes rainfall rate observations). Rice and Holmberg (Ref. 45) 
used a thunderstorm and nonthunderstorm characterization (and included 
geographic climatic regions) . 

It was necessary that meteorological quantities that might be 
used be available. Thus, while drop-size distribution aloft was the 
best quantity, it was unfortunately unavailable. 

The meteorological quantity most often used is R, the rain 
rate, which is both available and numerical. In addition to being re- 
lated to the reflectivity factor, Z, rain rate has another useful fea- 
ture - Information can be obtained on the probability of occurrence of 
various rain rates. It is this feature of R that makes it most useful 
in microwave communication, although this feature is of lesser interest 
to the radar meteorologist. There was, therefore, a third requirement 
for the meteorological data - that it be potentially capable of provid- 
ing probability estimates. 

For the purpose of microwave propagation and communication, 
there are specific expectations to be provided by a rain characteriza- 
tion. However, it is not a priori apparent that a specific characteri- 
zation of a rain type (say, nonthunderstorm) would optimally describe 
the statistics of the reflectivity factor, and the statistics of area, 
and the statistics of height distribution, and the statistics of inter- 
cell distance, etc. The question faced was if one characterization does 
not describe these well, what alternative characterization does, bearing 
in mind that only a practical scheme was acceptable. 

There are at least two approaches to a characterization. In 
one approach separate classes are identified such as thunderstorm and 
non thunder storm, or alternatively the classes: thunderstorm, warm 

front rain, cold front rain, squall line rain, etc., are identified and 
the statistics of each class determined. While the identification into 
classes did not meet the desire that it be numerical, it was neverthe- 
less acceptable. In the use of rain types or classes, it is obvious 
that the class definition be unambiguous. This was a real problem. 
Austin (Ref. 44) was careful to define what she meant by "thunderstorm 
rain." It had to meet the criterion that at least one of four observing 
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stations, within 40 mi, report thunder on the day that the rain occurred. 
In discussions sponsored by GCIR the term "excessive precipitation" is 
used to identify rain as "thunderstorm rain." There is a need for 
arriving at a common, acceptable definition. There are further compli- 
cations in the case where a rain meets the definition of two classes. ‘ 
Suppose there are some thunderstorms embedded in a squall line. Some' . 
priority would eventually be needed for assigning this rain to one 
particular class, of a decision reached on how to classify such cases. 
Another approach to icharacterization was to select a set of (plausible) 
meteorological parameters in addition to some of the above rain typing, 
and by regression or otherwise relate these to Z, if indeed there are 
significant correlations. The statistics would then be about some model 
relation between Z and other parameters. Probability estimates would 
use the probability of occurrence of different values of these parameters, 
singly or jointly as the model indicates. By this approach it may be 
possible to coalesce some rain types, e.g., thunderstorm and nonthunder- 
storm convective shower, into single types or the numerical equivalent 
of these types. 

In the statement that there are other "meteorological param- 
eters" of potential interest in categorizing rains, note that these were 
largely available from routine meteorology, and thus obtaining them did 
not pose a significant burden on the program. 

4.2 ELEMENTS IN RAIN CATEGORIZATION 

There was a body of knowledge available to aid in identifying 
rain classes and parameters. 'Local Influences which can impact on th'e 
analysis are also considered in this section. 

4.2.1 Physical Considerations 

The physics of precipitation, although incomplete, may “provide 
clues to the meteorological quantities more likely to be useful. • The 
phenomena that lead to the growth of raindrops are condensation on the 
precipitation particle and collision (accretion, coalescense) . In ad- 
dition, lifting of the moist air, which then cools, is needed. Upward 
air movements occur at different rates. In a warm front situation air 
is lifted slowly and drops grow principally by condensation. In a con- 
vective situation there are faster upward air movements and with drop- 
lets of different size and with different velocities there is an accen- 
tuated relative motion between the larger and smaller drops enhancing 
the growth of the larger drops by collision with the smaller drops. If, 
in addition, there is small scale turbulence the growth possibility by 
collision is further enhanced. There are processes that cause drops to 
decrease in size. There is evaporation and also drop break up if it 
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becomes so large that it distorts and breaks up, or if the outer layer 
freezes followed by the frozen shell breaking and releasing the drop 
interior as smaller drops. The probability of finding an ample supply 
of large drops, and thus higher Z, depends on having several growth 
processes occurring while drop diminution processes are of small impor- 
tance. 

Even though the foregoing discussion is simplistic it does 
point to categorization schemes, to the identification of desired param- 
eters and to potential concerns in interpretation of results. For aid 
in categorization there are data available from rawinsonde flights, from 
data transmitted on Weather Service circuits and from some special local 
measurements. The relationship between the requested information and 
.the physically significant information may be indirect; for example, an 
estimate of the relative intensity of the vertical velocity is obtained 
from the air "stability." There is a potential concern that the trans- 
formations between Z and R and attenuation and similar quantities pre- 
sume a certain drop size distribution. However, the discussion of the 
preceding paragraph suggests that the drop size distribution may differ, 
even possibly in the functional form, for different rain situations. 

Thus, we can intuitively expect more larger drops (higher Z's) in thunder- 
storms where there are pronounced vertical air flows, high moisture 
content in a cloud of large vertical extent and considerable small scale 
turbulence. In contrast, a hurricane may show the same rainfall rate 
but lower Z's as the small scale turbulence is less pronounced, and the 
repeated up and down motions that a drop undergoes in a thunderstorm is 
not common. Similarly, it may be possible to have the same Z's for 
different rainfall rates. 

The above discussion relating to rain characterization and the 
physics of precipitation are meant to point out some of the important 
features relating meteorological parameters with the statistical descrip- 
tion of precipitation as determined using radar. Although the radar 
analysis was to a large extent completed and documented in Volume II, 
a parallel effort was unfortunately not made in the meteorological area. 
Thus it yet remains to perform a definitive analysis to bring the meteor- 
ological aspects to a comparable position. 

4.2.2 Local Considerations 


In this section the meteorological situations which occur at 
Wallops Island and produce copious rains are identified. There are hur- 
ricanes and lesser storms of tropical origin; there are thunderstorms; 
there are coastal storms. The latter include lows that are near the 
coast, fronts that lie near the coastline and nearly parallel the coast; 
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such storms sometimes form as secondaries to a system located about 1000 
km to the west . (Other areas of the world have their own types of 
rainst — monsoon rains in India, coastal-orographic rains in Washington 
and Oregon.) In addition, there were, in the vicinity of Wallops Island, 
certain other local features. There is the ocean on the east, the 
relatively narrow (about 10 mi) strip of land to the west and the rela- 
tively broad Chesapeake Bay (about 30 ml) beyond. A convective type, 
storm. may need a continual thermal energy input from below to maintain 
itself. Passage over water could deny it an energy source. (There is 
still an energy source aloft from the heat release during condensation.) 
The. growth/decay history of a cell might depend on where the cell was 
observed and on its movement. The sea-breeze phenomena might have 
influenced the results as there was an added vertical velocity field. 
Another local consideration was the existence of a fairly heavily tra- 
velled airway nearby. The production of contrails (ice crystals) could 
have led to an intensification of rain by cloud seeding. 

4.2.3 Meteorological Descriptors 

a. Qualitative description . Individual days were described 
in a qualitative manner including general synoptic situa- 
tion. These descriptions contained such pertinent ele- 
ments as sea-breeze, presence of cirrus , and general 
storm movement. 

b. Rain rate . Quantitative measure of precipitation 
rate available from various rain gages. 

c. Vertical extent of cloud 

d. Horizontal extent of cloud 

- e. Stability measures . There are two stability measures 

currently ^distributed that are deducible from the radio- 
sonde flight, the lifted index and George's K-lndex. 

These were available from the local sonde, teletype 
reports FOUS 1 and 2 and FOUS 10 and 11, and the fac- 
; simile composite moisture chart. The K-index, for example, 

has been reported to give the highest correlation with 

^ radar activity, according to the National Weather Service 

Technical Procedures Bulletin, No. 71, March 6, 1972. 
i . Another class of stability measures is the vertical 

. velocity. This is not a measured quantity; it is com- 

puter-derived for the 700 mb level. It was obtained from 
the facsimile chart "Relative Humidity/Vertical Velocity," 
and from teletype reports FOUS 1 and 2, and FOUS 10 and 
'• 11 . 
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f. Storm motion . 

g. Moisture . Information on moisture-related quantitid's 
were available from several sources. There were the 
local and nearby sondes. Teletype reports FOUS 1 and 2 
gave the mean relative humidity of the lowest 3 levels of 
the 6 level primitive equation model (surface to about 
490 mb). It also gave the quantitative precipitation 
forecast. The composite moisture facsimile chart had ■ 
"precipitable water" and the average relative humidity ' 
(surface to 500 mb), each computed from radiosondes. The 
relative humidity /vertical velocity facsimile chart had ■ ’ ‘ 
the mean relative humidity as in FOUS 1. There was also ‘ 
the heavy precipitation forecast transmitted on the FOFAX 
circuit . 

h. Significant heights . There is a set of potentially sig- 
nificant heights/pressures from various sources. The use 
of pressures is motivated from thermodynamic considera- 
tions. The local and nearby radiosondes routinely report 
these. In addition the local sonde provided height and 
pressure of the -4°C level, the 0° level and the tropo- 
pause. The freezing level is also on the composite 
moisture facsimile chart, but it is noted that in the 
case of several freezing altitudes, only the lowest one 
is given. The tropopause is shown on the "Upper Levels/ 
Tropopause Data" facsimile chart. 

i. Winds . Winds from local measurements and from nearby 
stations are given to several levels. While it has been 
reported that in the mean, the echo drift velocity be- 
comes closest to the 10,000 ft wind, information for 
storms to varying heights may refine this. Consequently 
the analysts anticipated using the winds at 8,000, 10,000 
... 18,000 ft. It has also been reported that echoes 
move at a lower speed than the winds , with the lower 
speeds associated with the stronger echoes. This suggests 
an analysis which compares the difference between wind 
and echo drift with Z. Wind shear, the difference of 
wind between two altitudes, appears to have an effect on 
the vertical extent of shower clouds and storm intensity. 

j . Temperatures . With the thought that the higher the 
temperature the greater is the amount of water the ai:^ 
can hold and with the observation that as one moves 
equatorward, toward the warmer temperatures, the thun- 
derstorms extend to higher altitudes, (significantly. 
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there is a higher tropopause, too), the temperatures at 
selected standard levels are appropriate to the present 
experiment. These data come from radiosonde data and are 
given on standard level facsimile charts. A closely 
related quantity is the thickness of the surface to 500 
mb region. 

1 

5. SUMMARY OF DATA 

\ 

\ 

In this section are listed the data available from the various 
sources taken during .the 1972 and 1973 experiments. The radar data and 
that appearing in photographic form are -tabulated in Section 5.1. In 
the, meteorological aspects all the data sources and the data that have 
been cataloged and are available at APL for analysis are listed. These 
are found in Section 5.2. 

5.1 RADAR AND PHOTOGRAPHIC DATA 

This section lists in tabular form all the data available from 
the spatial statistics 1972-1973 experiment. 

For each day during which precipitation was observed the total 

niimber of each of the following is listed in Table 5.1. 

i 

1. Polaroid photographs. 

2. 360 deg scans. Also, the number of such scans on 
Polaroid (P) and on digital tape (T) . 

3. 2 deg elevation PPI scans. (P) and (T) as above. 

4. Raster scans. 

5. Sweeps in all PPI raster scans on that day. 

6. Average number of sweeps in each raster scan. 

7. Individual RHI scans. (P) and (T) as above. 

8. RHI raster scans. 

9. Sweeps in all RHI raster scans on that day. 

10. Searchlight mode observations. 
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Table 5.1 

Radar Data Available from the Spatial Statistics 
1972—1973 Experiment 


( 1 ) ( 2 ) 


Data 

ftihroid 

2 

Photos 

360 

P 

nns 

39 

15 

9/18 

95 

23 

9/19 

62 

17 

9/21 

26 

24 

9/25 

10 

10 

9/27 

154 

33 

9/29 

76 

30 

9/30 

76 

IS 

10/5 

9 

9 

10/6 

16 

6 

Total 



10 Days 

563 

182 

1973 


■1 

5/9 

5/14 

47 

62 

H 

1 1 

39 

84 

16 

14 

5/23 

71 

15 

5/24 

133 

17 

5/29 

184 

16 

6/13 

197 

17 

6/18 

196 

14 

6/20 

34 

15 

6/21 

209 

14 

6/22 

307 

19 

6/26 

22 

10 

6/27 

113 

16 

6/28 

45 

15 

6/29 

353 

16 

7/3 

198 

16 

7/5 

55 

12 

7/10 

281 

13 

7/11 

175 

11 

7/18 

33 

12 

7/20 

60 

14 

7/27 

163 

14 

7/31 

116 

13 

8/1 

225 

12 

8/2 

171 

14 

8/3 

225 

17 

8/13 

277 

13 

8/14 

135 

17 

8/15 

200 

16 

8/16 

147 

15 

8/20 

109 

15 

8/21 

154 

16 

Total 
33 Days 

4822 

480 



P = Polaroid 
T«Tape 
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In addition, Table 5.1 contains: 

11. Digital tape identification numbers. 

12. Video tape identification numbers. 

13. PPI photograph identification numbers recorded on 35 mm 
film, j 

Totals of all rain days in 1972 and 1973 appear on the last 
line for each year. 

In an earlier discussion it was pointed out that digital tape 
equipment was not available at the start of the experiment. Thus, under 
1972, no identification numbers appear for digital tapes. From 9 May to 
24 May 1973 are listed DV Identification numbers which mean digital 
tapes made from video tapes in "off-line" mode. Tape identification not 
bearing a D or DV refer to video tapes. 

5.1.1 Radar Data Analysis (Volume II) 

j 

This report (Volume I) on the rain statistics experiment has 
dealt with various aspects of the experimental program and does not 
contain any significant discussion on data analysis. However, the an- 
alysis program was an important but parallel effort which had serious 
Impact on the form of the experiment. Development of the concepts of 
the analysis program may be found in Ref. 5, 46, 47, 48 and 49. 

In the course of performing the analysis, as is the case in 
most such programs, changes were made to the program as it progressed. 
Data reduction and analysis techniques and the statistical results as 
finally evolved can be found in Volume II of this report. 

5.2 METEOROLOGICAL DATA 

5.2.1 Sources of Data 


The meteorological data have come from several sources: 

a. Data generated by the Weather Service Office (WSO) at 
Wallops Station. 

b. Teletype Service A (Flight Services) circuit 005 reports 
were provided by Wallops WSO. 

c. Telet)^e Service C (Weather Service) circuit 31 reports 
were provided by Wallops WSO. 
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d. A teletype RAWARC circuit report was provided by Wallops 
WSO. 

e. Facsimile Charts (National Schedule) were provided by 
Wallops WSO. 

f. Fascimile Charts (Forecast Schedule) were obtained from 
the National Meteorological Center, Sultland, Maryland. 

g. Fascimile Weather Radar maps transmitted by land lines 
from Patuxent River, Maryland were received at SPANDAR. 

h. Radar weather forms prepared at Patuxent River, Maryland. 

i. Data generated by APL personnel. 

J . Hourly precipitation data for Maryland and Virginia lo- 
cations issued monthly by the U.S. Department of Com- 
merce, Asheville, North Carolina. 

k. Climatological data for Maryland and Virginia locations, 
issued monthly by the U.S. Department of Commerce, Ashe- 
ville, North Carolina. 

l. Daily Weather Maps, Weekly Series, issued by the U.S. 
Department of Commerce and prepared at Suitland, Mary- 
land. 

m. Ocean Water Surface Temperature charts, issued monthly 
by the U.S. Coast Guard (Department of Transportation). 

5.2.2 Data Catalog 

Data have been extracted from the above sources and have been 
filed and are readily retrievable. The data are listed next; the source 
is identified parenthetically using identifiers in 5.2.1. (The numbers 

5.2 in the listing of the above paragraph will be omitted for simpli- 
fication.) 

a. Wallops Island Rawinsonde evaluations of OOOOZ and 1200Z 
soundings (l.a). 

(1) Point values of height, pressure, temperature and 
dew-point depression reported at significant and 
mandatory levels. 

(2) Showalter Stability index. 
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(3) Vertically averaged wind speed and direction re- 
ported at approximately 1000 f t^intervals. 

b. Wallops Island OOOOZ rawinsonde evaluations at 1000 ft 
levels for periods of SPANDAR operations (l.a). 

(1) Temperature, index of refraction, pressure, den- 
sity, dew-point, relative humidity, absolute humid- 
ity, wind speed and direction. 

c. Nearby rawinsonde data from stations within about 1000 km 
of Wallops Island (l.c). These contain same information 
as 2. a, above. Data for lower altitudes (under 50,000 
ft) were filed. 

d. Wallops Island Surface Observations form lOA (l.a). 

These give state of sky, weather, celling, visibility, 
sea level pressure, temperature, dew-point and surface 
wind speed and direction. These are given hourly from 
1200Z to 2100Z, but; on weekends are given at 1500Z, 

1800Z, and 2100Z. 

e. Wallops Island Surface Observations form lOB (l.a). 

This form supplements form lOA and gives time of begin- 
ning and ending of rain, showers, thunder, etc., and 
hourly amounts of precipitation. 

f. Wallops Island rain gage charts (l.a). Weighing gage 
with time resolution about 3-4 min and rain resolution 
of .01 or .02 in. 

g. Charts from rain gage located at Snow Hill VOR (1.1). 
Tipping bucket gage, each tip indicates 0.2 mm rain, 
with time resolution of about 0.1 to 0.2 min. Total 
rainfall measured separately. 

h. Surface weather data from about 500 km radius (l.b). 

These are observations taken at airport stations. The 
data that were filed covered (1) periods when stations 
within 140 km of SPANDAR reported rain, and/or (2) 
periods when SPANDAR was taking data on rain. These 
data include observations of rain including intensity 
and thunder. There are also other meteorological sur- 
face observations of rain including intensity and 
thunder. There are also other meteorological surface 
observations in the airport reports. 


66 



1 . 


Surface weather maps issued at 3-hour intervals (l.e). 

Daily surface weather maps (1.1). 

j. Wdather map analyses and summaries, particularly for 
rain periods (l.b and l.c). 

k. Upper Air Maps (l.e). These maps for the 850, 700, 500, 

300, and 200 mb levels Include wind, temperature, dew- 
point depression and heights and are issued twice daily. 

Also separate wind maps. 

l. Weather Depiction charts (l.e.). These maps, issued 
eight times daily, show sky cover and ceiling. 

m. Patuxent Radar facsimile charts (l.g). These transmis- 

sions were recorded three times per hour during the work- 
ing day and show the PPI scope with annotations. There \ 

are also the radar weather forms (l.h). 

n. Radar Depiction Charts (l.e). These charts. Issued 14 
times dally, give the combined analysis from the network 
of radar stations. Echo drift velocities are given; 
rain categorization and intensities are given. 

o. Coded radar reports (l.b). These reports give the areas 
of returns with fractions of area covered, separate cells 
and their diameters. There is also an estimate of rain 
Intensity, drift velocity and maximum heights. Echo 
intensities are given. The data retained were those that 
indicated rain in the SPANDAR area. Both Patuxent and 
Atlantic City radar reports are in the list of reporting 
stations. 

p. Composite Moisture index chart (l.e). Issued twice 

daily, this contains four charts: a stability chart 

which contains both "lif ted-lndex" and "K-lndex"; pre- 
cipltable water (surface to 500 mb) ; freezing level 
(lowest level that 0°C is reported if there is more than 
one level) ; and the average relative humidity (surface to 
500 mb arid weighted by pressure depth) . 

q. Relative Humidity /Vertical Velocity (l.e). Issued twice 
dally. The dynamically derived vertical velocities are the 
700 mb (about 10,000 ft) values. 
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r. Lifted Index chart (l.f). 

s. Surface and 1000 to 500 mb thickness chart (l.e). 

t. Vortlcity charts (l.e). Both barotroplc and baroclinic. 

u. Three-dimensional trajectory forecasts (l.c). For 
selected cities. Also has K-index. • 

V. Climatological data (l.k). This Includes the 24-hr 

precipitation as measured by cooperative observers as 
• well as by the regular weather stations. There are 

about ten locations within 40 mi of SPANDAR and about 
double that within the SPANDAR observing range. 

w. Hourly precipitation data (l.j). These contain reports 
from six stations with recording rain gages within 
the SPANDAR observing range. . . 

X. Water temperatures (l.d and l.m). 

y. Miscellaneous Weather Service outputs (l.b,c,e,f). 

Synoptic discussions. Numerical products. 

5.2.3 Data Reduction 

Detailed tabulation and data reduction of the pertinent meteor- 
ological data are in Refs. 50, 51, 52, 53, arid 54/' These are considered 
as part of this Final Report and are included in Volume III as Appendix 
D.. 


5.2.4 1973 Rain Comparison . . - 

The total monthly rainfall in the summer of 1973 was compared 
with the 10-year period 1962-1971, Ref. 55.’ Table 5.2 gives for each of 
the summer months, the 10-year average, the 1973 average, the wettest 
and driest month and climatology based on an even longer period but 
using fewer stations. ' 

Table 5.3 gives the average number of days of different rain- 
fall amounts, based on the same 10-year period compared with that of 
1973. . , 

The tables show that the decision to have the 4 -months of 
observation extend from May to August was. a good one. The month of May 
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Table 5.2 


Monthly Rainfall (in.) 



May 

Jun 

Jul 

Aug 

Sep 

10-year average 

mm 

■a 

3.8 

4.1 

3.3 

1973 average 

3.7 

3.7 

3.2 

■B 

2.5 

Average, wettest month 

msM 

5.6 

(1963) 

6.8 

(1969) 

8.7 

(1967) 

H19 

Average, driest month 

0.6 

(1965) 

—■ 


mi 


Climatology 

3.6 

3.6 

■B 

5.6 

4.0 


Table 5.3 


Average Number of Days of Specified Rainfall 


Rainfall Equals 
or Exceeds 
(in.) 

May 

Jun 

Jul 

Aug 

Sep 

Av 

1973 

Av 

1973 

Av 

1973 

Av 

1973 

Av 

1973 

.01 

8.3 

12.1 

m 

8.7 

m 

8.3 

7.8 

10.5 

6.0 


.1 

m 

8.7 

5.3 



5.5 

5.5 

8.2 

4.3 

ED 

.5 

1.8 

2.9 

2.3 

2.6 

2.6 

2.6 

2.3 

m 

2.2 

1.6 

1.0 

0.6 

0.4 

0.9 

0.9 

1.0 

0.7 

1.2 

2.0 

1.0 

0.9 


was rainier than the 10-year average and about the same as the longer 
period climatological value. There were substantially more days of 
measurable rain. 

The month of June was about average with respect to amount of 
rain and number of rainy days. The month of July was drier than normal 
in amount of rain and somewhat below normal in terms of days of rain. 
August had substantially more rain and more days of rain. In total, the 
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sunnner observing months of 1973 had more rain than the 10-year average 
indicated and more days of rain, but it was no record-breaker. In terms 
of the total rainfall for the 10-year period and the 4 months of obser- 
vation — a total of 40 months — 1973 was wetter 60% of the time, drier 32% 
of the time and about the same (within 0.25 in. of rain) 8 % of the 
time. 
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